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SECTION  I 


INTRODUCTION 

Bounded-wave  electromagnetic-pulse  (EMP)  simulators  utilize  a  pair  of 
conductors  to  guide  and  confine  the  electromagnetic  fields  in  the  test 
volume.  A  portion  of  the  conductors  is  often  in  the  form  of  two  parallel, 
rectangular  plates  or  strips.  Figure  1  shows  a  typical  geometry  of  the 
bounded-wave  simulator.  It  consists  of  two  identical,  long  conductors 
joined  together  at  one  end  to  an  electromagnetic  pulser,  and  at  the  other 
end  to  a  load  termination.  This  conductor  configuration  can  be  regarded 
as  made  up  of  a  section  of  a  finite-width,  parallel-plate  transmission  line 
sandwiched  between  two  conical  transmission  lines.  The  test  volume  is 
located  between  the  parallel  plates. 

When  the  pulser  is  fired,  electric  currents  rush  out  from  it  onto 
the  two  conductors.  These  currents  are  accompanied  by  electromagnetic 
fields  which  propagate  from  the  pulser  toward  the  load.  Because  of  the 
special  geometrical  design  of  the  conductors,  the  electromagnetic  fields 
will  be  predominantly  in  the  transverse-electromagnetic  (TEM)  mode.  This 
feature  is  highly  desirable  since  the  TEM  mode  closely  resembles  the  plane- 
wave  nature  of  the  actual  nuclear  EMP.  The  load  is  suitably  chosen  to 
minimize  reflection  from  the  termination. 

One  approach  to  achieving  an  overall  understanding  of  the  performance 
of  the  bounded-wave  simulator  consists  in  isolating  and  analyzing  specific 
aspects  of  the  simulator  operation.  Such  aspects  include  the  electromagnetic 
propagation  characteristics  of  the  parallel-plate  and  conical  transmission 
lines.  These  have  been  studied  in  references  1,  2  and  3.  Another  aspect 
of  special  importance  is  the  propagation  of  the  electromagnetic  excitation 
from  the  conical  transmission  line  onto  the  parallel  plates.  The  manner 
of  transition  of  the  excitation  across  the  junction  directly  affects  the 
electromagnetic  fields  in  the  test  volume  between  the  parallel  plates. 

This  report  presents  a  theoretical  study  of  the  simulator  junction. 

The  analysis  is  based  on  a  geometrical  model  shown  in  figure  2.  The  model 
consists  of  a  semi-infinite  section  of  a  parallel-plate  transmission  line 
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odel  for  Studying  the  Junction  in  a  Bounded-Wave  EMP  Simulator. 


connected  to  a  finite  section  of  a  conical  transmission  line.  Both  trans¬ 
mission  lines  are  assumed  to  be  made  out  of  infinitely  thin  and  perfectly 
conducting  sheets.  The  conical  transmission  line  is  excited  at  the  apex 
by  a  given  pulser.  It  is  required  to  calculate  the  resulting  electromagnetic 
excitation  propagating  onto  the  parallel-plate  section. 

The  model  clearly  reproduces  exactly  the  geometry  of  the  pulser  end 
of  the  simulator.  One  has,  however,  exploited  the  weakness  of  the  reflection 
from  the  termination  to  replace  the  load  end  with  a  semi-infinite  continuation 
of  the  parallel-plate  section.  The  geometry  becomes  sufficiently  simple  to 
permit  an  approximate  analytical  solution. 
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SECTION  II 


GEOMETRY 


In  studying  the  bounded-wave  simulator  model  of  figure  2,  it  is 
mathematically  convenient  to  consider  instead  an  equivalent  model  consisting 
of  a  single  conductor  over  an  infinite,  conducting  ground  plane.  The  image 
of  the  conductor  below  the  ground  plane  replaces  the  second  conductor.  The 
apex  of  the  conductor  in  this  equivalent  model  is  to  be  driven  against  the 
ground  plane  by  an  equivalent  pulser. 

One  sets  up  a  rectangular  coordinate  system  such  that  the  x-y  plane 
coincides  with  the  conducting  ground  plane,  as  shown  in  figure  3.  Let  the 
conductor  be  situated  in  the  upper  half-space  z  >  0.  Its  apex  is  located 
at  the  coordinate  origin.  Its  semi-infinite-strip  section  is  parallel  to 
the  ground  plane  along  the  positive  x-axis.  Let  the  width  of  the  strip 
be  denoted  by  2a  and  its  height  above  the  ground  plane  by  h.  The  conical 
section  is  actually  an  isosceles  triangle  of  height  d  and  base  length  2a. 

It  is  inclined  to  the  ground  plane  at  an  angle  0  such  that 

sin  0  =  j  (1) 

It  will  be  convenient  to  introduce  at  the  same  time  a  two- 
dimensional  rectangular  coordinate  system  on  the  surface  of  the  conductor 
itself.  In  this  system  a  point  on  the  conductor  is  uniquely  defined  by 
a  pair  of  coordinates  (£,y),  as  shown  in  figure  4.  Here  y  is  the  same  as 
the  y-coordinate  of  the  three-dimensional  rectangular  coordinate  system, 
while  S  is  related  to  x  and  z  by  linear  relations.  Let  the  triangular 
plate  be  denoted  by  and  the  strip  by  S2*  The  entire  surface  of  the 
conductor  will  be  denoted  by  S  (=  S^+S^)*  Then  these  linear  relations  are: 

1.  On  S1  (0  <  5  <  d ) 

x  *  £  cos  9  (2) 

z  ■  5  sin  9 
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Figure  3.  Three-Dimensional  Coordinate  System  for  the  Bounded-Wave 
Simulator  Model. 
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Figure  4.  Two-Dimensional  Coordinate  System  for  the  Surface  of  the  Bounded-Wave  Simulator  Model. 


2.  On  S2  (d  <  C  <  “) 

x  ■  £-d  +  d  cos  9  (3) 

z  -  h 

One  may  call  £  the  longitudinal  coordinate  and  y  the  transverse  coordinate. 
Note  that  the  apex  of  the  conductor  is  situated  at  £  =  y  *  0. 
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SECTION  III 


ELECTROMAGNETIC  POTENTIALS 


The  bounded-wave  simulator  model  of  figure  3  can  be  analyzed  In  the 
frequency  domain.  One  considers  the  pulser  as  a  time-harmonic  voltage 
source  driving  the  apex  of  the  conductor  against  the  ground.  Its  variation 
with  time  t  is  described  by  the  time  factor  e^ut  with  angular  frequency  u. 

For  simplicity  it  will  be  assumed  that  w  is  positive.  Results  corresponding 
to  negative  u>  can  be  obtained  from  those  for  positive  u  by  complex  conjugation. 
Under  this  excitation  a  time-harmonic  electric  current  flows  on  the  conductor 
as  well  as  on  the  ground.  It  suffices  for  the  analysis  to  calculate  the 
surface  current  density  on  the  conductor. 

1 

In  this  study  the  current  on  the  conductor  will  be  calculated  by 
using  the  electric-field  integrodifferential  equation.  This  approach  has 
great  methodological  appeal  in  that  it  enables  one  to  obtain  an  exact 
formulation  of  the  problem,  despite  the  difficult  geometry  of  the  conductor. 
This  exact  formulation  serves  as  a  solid  base  from  which  approximations  can 
safely  proceed.  The  idea  of  the  formulation  is  to  relate  the  electric 
field  outside  the  conductor  to  the  current  on  the  conductor,  and  then  impose 
boundary  conditions  on  the  electric  field  at  the  surface  of  the  conductor. 

For  this  purpose  it  is  convenient  to  introduce  the  electromagnetic  potentials. 

Let  £  denote  a  point  in  the  upper  half-space  z  >  0: 

r  =■  (x,y,z)  (4) 

and  let  _E(ir)  and  B(r)  be  the  electromagnetic  fields  at  x_  generated  by  the 
currents  and  charges  on  the  conductor  and  the  ground  plane.  Then  E  and  B 
can  be  derived  from  a  pair  of  scalar  and  vector  electromagnetic  potentials 
9 (r_)  and  A(r) : 


E(r)  -  -Vcp(r)  -  jwA(r) 
B(r)  -  V  x  A(r) 


(5) 
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In  Che  Lorentz  gauge  these  potentials  individually  satisfy  the  Helmholtz 
equation: 


(v2  +  lc2)<p(r)  -  0 

(72  +  k2  / A (r)  -  0 

where 

ic  - 

o  o 


(6) 


(7) 


and  cq  and  uQ  are  the  permittivity  and  permeability  of  free  space.  By  use 
of  suitable  Green's  functions  the  solutions  of  equations  6  can  be  represented 
as  integrals  over  the  surface  charge  density  a  and  the  surface  current  density 
K  on  the  surface  S  of  the  conductor  alone: 


9(£) 


G^r.r'Wr'JdS' 

V£’£,)Kx<£,)ds' 


Ay(r) 

Az(r) 


S 


Gj,(£,»£.,)Ky  (£,'  )ds ' 


V^')Kz(£')dS’ 


(8) 


In  these  expressions  jr'  denotes  a  point  on  the  surface  S: 


r'  -  (x',y’,z’) 


(9) 


and  dS'  is  a  surface  element  of  S  at  _r'  . 

The  functions  G_  and  G„  are  Green's  functions  of  the  Helmholtz 
D 

equation  in  the  half-space  z  >  0: 


They  satisfy,  respectively,  the  homogeneous  Dirichlet  and  Neumann  boundary 
conditions  on  the  ground  plane  z  «  0 : 


GjjOr ,r ’ )  -  0 


IF  *  0 


(ii) 


By  introducing  the  image  r"  of  r_'  with  respect  to  the  ground  plane,  such 
that 


r"  -  (x'.y'.-z*) 


(12) 


one  can  construct  G_  and  G„  from  the  free-space  Green's  function  G  : 

D  N  r  o 


The  results  are: 


Go^’> 


e-jk|r-r’  | 
4ir|r  -  r'  | 


(13) 


Gd (r,r')  -  Go(r,r’>  -  GQ(r,r") 
GN(r ,r ’ )  -  Go(r,r')  +  Go(r,r") 


(14) 


The  free-space  Green's  function  Gq  in  equation  13  already  has  a  simple, 
closed  form.  Nevertheless,  in  later  calculations,  it  will  be  advantageous 
to  work  also  with  its  Fourier- integral  representation: 

-j Io(x- x')  +  S(y -y')  +  y(z  -  z')] 

G  (r,r') - -  da  dB  dy  - - 2 - 5 - 5 - 2 -  (15) 

0  (2*r  J-w  a  +B  +y  -k 

For  this  expression  to  be  identical  to  equation  13,  it  must  be  understood 
that  k  has  an  infinitesimal,  negative  imaginary  part: 

k  -*■  k  -  je  e>0  (16) 
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The  surface  integrals  in  equations  8  take  on  a  simpler  form  if  one  goes 
over  from  the  three-dimensional  coordinate  system  (x,y,z)  to  the  two- 
dimensional  coordinate  system  (£,y)  set  up  on  the  conductor's  surface  S, 
as  defined  in  equations  2  and  3.  In  the  latter  coordinate  system  the 
surface  current  density  K  has  only  two  components  and  Ky.  The  components 
Kx  and  in  the  former  coordinate  system  are  both  related  to  K^. 


1.  On 

K  -  K.cos  9  (17) 

x  5 

K  -  K.sin  6 

z  e 

2 .  On  S2 

-  Ee  (18) 


In  terms  of  and  K^,  the  vector  potential  A  reads 


A  (r)  -  u 
x  —  o 


GD(r,r')cos  9  K  (r')dS'  +  uo  f  GD(r,r')K  (r')dS' 


Ay<I>  -  ^ 


GD(i’^,)Ky(i')dS' 


(19) 


a2(£)  *  Uo  (  GN^-,-^sin  9  ^(r')ds' 


In  these  expressions  the  components  of  it'  are  related  to  the  source  coordi¬ 
nates  (C' ,y ' )  according  to  equations  2  and  3,  while 


dS’  -  dC'dy' 


(20) 


18 


I 

I 

I 
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SECTION  IV 


EXACT  FORMULATION 

There  are  two  electric-field  integrodifferential  equations  for  this 
problem.  They  are  formulated  by  imposing  appropriate  boundary  conditions 
on  the  two  tangential  components  E, (r)  and  E^(r)  of  the  electric  field  at 
every  point  r  on  the  conductor's  surface  S.  To  derive  these  boundary 
conditions  one  regards  the  pulser  as  the  supplier  and  maintainer  of  a 
definite  tangential  electric  field  with  components  E°(r)  and  E°(r)  on  the 
conductor.  One  then  requires  that  E^(r)  and  E^O:)  assume  these  prescribed 
values  on  S,  that  is. 


E  (r)  =  E°(r) 
Ey(r)  =  E°(r) 


(23) 


for  _r  lying  on  S. 

In  its  idealized  form  the  pulser  is  a  localized  voltage  source.  It 
generates  a  finite,  time-harmonic  voltage  VQ  across  the  infinitesimal  gap 
between  the  two  triangular  plates  at  the  apex  of  the  conical  transmission 
line  in  figure  2.  The  tangential  components  of  the  pulser  electric  field 
are  therefore  given  by 


E°(r)  -  -  Vo5(5) 

E°(r)  =  0 

y  - 

The  5-function  is  normalized  to  make 

J  d£  5(5)  “  1 


(24) 


(25) 


Therefore  the  voltage  across  the  apex  of  the  conductor  and  the  ground  plane 
in  figure  4  is  y  Vq  since 
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'  1  “5  E°(ri  -  f  V0 


The  integration  over  £  from  0  to  00  picks  up  only  one  half  of  the  contribution 
from  the  5-function  in  E°. 

Combining  equations  5,  23  and  24  one  obtains  the  following  pair  of 
boundary  conditions  for  the  potentials : 


Je  cp(r)  +  jui  A^(r) 


V<« 


+  *,(£> 


for  r  lying  on  S.  The  component  A  is  related  to  A  and  A  as  follows: 

^  X  Z 


1.  On  S, 


A^  =  Axcos  0  +  Azsin  0 


2 .  On 


A_  ”  A 
S  x 


Substituting  equations  28  and  29  into  equations  27  and  eliminating  the 
potentials  by  use  of  equations  19  and  22,  one  transforms  the  boundary 
conditions  into  a  pair  of  coupled  integrodif ferential  equations  for  the 
two  components  and  of  the  surface  current  density: 

fiF  GpOr.r^V'-KO^dS'  +  k2  j  GH(r  ,r ’  )K,  (r '  )dS' 


-  J<»eoVo5(0 


f-  GjjCr.r'jV'-KCr'JdS’  +  k2  f  GD(r,r’)K  (r’)dS’  -  0 
y  'c  * 
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The  quantity  Gg(r_,r/ )  is  a  hybrid  Green's  function.  It  is  a  linear 

combination  of  the  Dirichlet  and  Neumann  Green's  functions  G_,  and  G„,  and 

D  N 

is  defined  in  the  following  manner: 


1.  For  r  and  r’  on  S. 


GH **  GD(r..r')cos29  +  GN(r,r.')sin28 


(31) 


2.  For  r_  on  and  r_'  on  or  vice  versa 


Gjj(r^,r_’ )  «  GD(r,r')cos  6 


(32) 


3.  For  r  and  r*  on  S, 


*  Va-L"> 


(33) 


Equations  30  are  in  part  differential  equations.  They  must  be  supple¬ 
mented  by  additional  boundary  conditions.  These  can  be  taken  as  the 
requirement  that  the  components  of  the  surface  current  density  normal  to 
the  edges  of  the  conductor  be  zero.  On  the  edges  are  given  by  y  =  ±(a/d)£. 
On  S2  they  are  given  by  y ■ ±  a.  The  boundary  conditions  on  these  edges  are: 


1.  On 


K  *  ±  f  K_ 
y  d  5 


for  y  =  ±  -j-  5 


(34) 


2.  On  S2 

K  =  0  for  y  ■  ±  a 

y  3 


(35) 


The  integrodif ferential  equations  30  and  the  boundary  conditions  34  and  35 
provide  an  exact  formulation  of  the  bounded-wave  simulator  model  of  figure  2. 


SECTION  V 


APPROXIMATE  FORMULATION 


Equations  30  are  a  pair  of  coupled,  two-dimensional,  partial  integro- 
differential  equations.  Their  solution  constitutes  a  formidable  mathematical 
problem.  At  the  moment  it  does  not  appear  possible  to  forge  ahead  with  the 
solution  without  introducing  approximations.  In  the  following  calculations 
a  number  of  approximations  will  be  made.  They  enable  one  to  obtain  an 
approximate  solution  of  the  problem  in  analytical  form. 

One  type  of  approximation  to  be  applied  is  aimed  at  simplifying  the 
basic  equations  of  the  problem.  What  one  essentially  does  is  to  smooth  out 
the  transverse  variations  of  the  surface  current  density  by  an  averaging 
process.  In  this  way  one  is  led  to  an  approximate  integrodif ferential 
equation  for  the  total  longitudinal  current  flowing  on  the  conductor. 

This  equation  constitutes  an  approximate  formulation  of  the  bounded-wave 
simulator  problem. 

Consider  the  first  member  of  the  two  basic  equations  30 


—  [X(C  )  dy’GD(c,y,e\y')  [f^T^U’.y’)  +  f^.  Ky(S' ,y')  j 


3C 


+  k 


d£’ 


•A  (£’> 


dy'GgU.y.S’.yMKgtt’.y')  =  -  (36) 


The  quantity  A(£')  appearing  in  the  limits  of  the  y '-integration  is  defined 
as  follows: 


HV) 


aC'/d 


0  <  V  <  d 
d  <  £'  <  oo 


(37) 


The  difference  in  the  functional  form  of  \(£')  for  the  two  ranges  of  £' 
reflects  the  fact  that  for  £'  <  d  the  conductor  is  a  triangular  plate,  while 
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-  Tgg'- 


for  V  >  d  it  is  a  rectangular  strip.  The  left-hand  side  of  equation  36 
is  proportional  to  the  5-component  of  the  electric  field  on  the  conductor 
generated  by  the  surface  current  and  charge  densities  on  the  conductor. 
The  kernels  Gp  and  Gy  can  be  regarded  as  coupling  functions  relating  the 
surface  current  and  charge  densities  at  the  source  point  (C',y')  to  the 
longitudinal  component  of  the  electric  field  at  the  field  point  (5,y). 

One  can  extract  from  equation  36  an  approximate  equation  for  the  total 
longitudinal  current  by  first  neglecting  the  transverse  surface  current 
density  and  then  averaging  out  the  dependences  in  the  two  kernels  Gy 
and  Gy  on  the  transverse  coordinates  y  and  y'. 

The  two  specific  approximations  can  be  justified  on  the  ground  that 
one  has  under  consideration  a  very  long  conductor  which  is  driven  at  one 
end  by  a  voltage  source.  In  this  situation  the  surface  current  will  be 
flowing  predominantly  in  the  longitudinal  direction.  One  can  expect  the 
following  inequality  to  hold  generally: 


\\\  >>  |Kyl 


The  discarding  of  iC^  from  equation  36  will  therefore  incur  only  a  small 
error.  Furthermore,  precisely  because  the  conductor  is  long,  it  will  be 
a  good  approximation  to  replace  the  two  coupling  functions  G_  and  Gu  by 
their  averages  in  the  transverse  direction.  That  is  to  say,  one  lets 


GD(C,y,f',y') 


GH(S,y,5',y') 


rD(C»£')  =  4x(S)X(C') 


MO  MO) 

j  dy  < 

-*(5)  -HO) 


dy,GD(5,y,C,,y') 


=  4x(s)A(C')  {^^dy  j  f  ^dy'GH(5'y*?' ,y,) 


HO  rU5') 


In  this  way  equation  36  becomes  an  ordinary  integrodif ferential  equation 
for  the  total  longitudinal  current  1(5): 


d_ 


deTD(5,5')  frr  I(C')  +  *2  f 


dc’rH(5,5')i(C) 


-jWeoVo6(0 


(41) 


where 


tHi) 

1(0  =  J  dy  K  (5,y) 
-A(0 


(42) 


Equation  41  is  a  closed  equation  for  the  total  longitudinal  current. 
Although  by  no  means  trivial,  it  is  significantly  simpler  than  the  pair  of 
exact  equations  30.  The  calculations  in  the  following  sections  will  be 
based  on  this  approximate  formulation  of  the  bounded-wave  simulator  problem. 
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SECTION  VI 


REDUCTION  TO  INTEGRAL  EQUATION 

The  integrodif ferential  equation  41  can  be  reduced  to  a  pure  integral 
equation  by  Fourier  transformation.  Introduce  the  double  Fourier  sine  and 
cosine  transforms  FQ  and  FH  of  the  kernels  Tq  and 


FD(p,p')  =*  P  dc  r  d^'rD(C,^')sin  pC  sin  p'C' 

v  ;  Jo  j0 


Vp’p,)  =(f)2 


•  OP  fCO 

dC 

0 


dC ’ rH(C » €  * )cos  pC  cos  p'C' 


(43) 


The  inverse  transforms  are 


rD(C,§') 


dp  |  dp’FD(p,p')sin  pC  sin  p'C’ 


(44) 


rH(c,e')  = 


dp 


dp'FH(p,p')cos  pC  cos  p'C' 


Substituting  equations  44  into  equation  41,  one  obtains 


d 

J, 


fCO  foa 

dq 


dp'Fn(q,p' )sin  qC  P(p')  +  k  dq 
J0  0 


dp’F  (q.p’)cos  qC  Q(p’) 


where 


=  -jmeoVo6(C) 


,00 


P(p') 


dC'sin  p'C'  1(5') 


(45) 


(46) 


Q(p')  -  f  dC 'cos  p’C’  I(C') 
Jo 


For  these  two  integrals  to  exist  one  must  assume  that  the  conductor  is  slightly 
resistive,  so  that  the  total  current  I  and  its  derivatives  vanish  at  infinity. 
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Next,  multiply  both  sides  of  equation  45  by  cos  pC  and  integrate 
over  C  from  0  to  00 .  The  operator  d/dC  can  be  handled  by  integration  by 
parts.  Then,  making  use  of  the  identities 


d£  sin  pC  sin  p'C  =  6(p-p') 


dC  cos  p £  cos  p'£  =  6(p-p') 
valid  for  p,  p'  >  0,  one  obtains 


(4; 


dp'  ^ pFD(p,p’)P(p')  +  k2FH(p,p')Q(p') J  = 


jwe  V 
oo 


(48 


The  two  functions  P  and  Q  are  closely  related  to  the  Fourier  cosine  trans¬ 
form  A(p)  of  the  total  current  I: 


A(p)  *  — 

TT 


d£  1(5) cos  p£ 


with  the  inverse  transform 


1(C)  =  dp  A(p)cos  pC 


From  equations  46  one  easily  sees  that 


P(p') 


-  f  p'A(p') 


Q(p')  =  j  A(p') 


Therefore,  substituting  equations  51  into  equation  48,  one  obtains 


I"  dp'  [-pp'FD(p,p’)  +  k2FH(p,p')j  A(p') 


2j(DE  V 
oo 


(49 


(50 


(51 


(52 


2: 


Equation  52  has  the  appearance  of  a  Fredholm  integral  equation  of  the 
first  kind  for  the  unknown  function  A(p').  Actually  it  is  an  integral 
equation  of  the  second  kind.  In  the  following  section  it  will  be  shown 
that  both  Fn  and  Fu  in  the  kernel  contain  a  6-function: 


FD(p,p')  *  P(p)5(p-p')  +  fD(p,p’) 
fh(P.P')  =  p(p)5(p-p')  +  fH(p,p’) 


(53) 


The  coefficient  of  the  6-function  is  the  same  in  both  cases.  Upon  substitution 
of  equations  53  into  equation  52,  the  latter  becomes 

j*oo 

(k2  -  p2)p(p)A(p)  +  dp' A(p,p')A(p') 

,l0 

where  the  new  kernel  A(p,p')  is  given  by 

A(p.p')  «  k2fH(p,p*)  -  pp'fD(p,p')  (55) 


2jwe  V 
0  0 


(54) 


Equation  54  is  an  exact  consequence  of  equation  41.  Being  essentially 
an  integral  equation  of  the  second  kind,  it  is  a  convenient  starting  point 
for  constructing  an  approximate  analytical  solution. 
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SECTION  VII 


EVALUATION  OF  FOURIER  TRANSFORMS 


Before  proceeding  with  the  solution  of  equation  54,  it  is  necessary  to 

obtain  explicit  expressions  for  the  two  quantities  p(p)  and  A(p,p').  This 

is  done  bv  calculating  the  double  Fourier  transforms  F_(p,p')  and  F  (p,p') 

D  n 

defined  in  equations  43. 

The  evaluation  of  the  two  double  integrals  in  equations  43  is  complicated 

by  the  fact  that  the  kernels  r  (£,£')  and  r  (£,£')  have  different  functional 

D  n 

dependences  on  their  arguments  over  different  parts  of  the  domain  of 

integration.  This  situation  can  be  appreciated  by  examining  the  £-£' 

plane  in  figure  5.  The  double  integrals  are  over  the  first  quadrant  of  the 

plane.  This  quadrant  has  been  divided  into  four  areas  labeled  1,2,3  and  4, 

with  the  boundaries  located  at  5  * d  and  £'  =d.  Across  these  boundaries 

the  functional  relations  between  (x,z)  and  £,  and  between  (x',z')  and  £', 

undergo  a  change  according  to  equations  2  and  3.  The  change  produces  a 

corresponding  discontinuity  in  the  functional  forms  of  T  and  r  .  This 

D  rl 

discontinuity  is  of  course  a  reflection  of  the  geometrical  discontinuity 
at  the  junction  of  the  triangular  plate  and  the  rectangular  strip. 

The  kernels  r_(£,£')  and  ra(£,£')  can  be  regarded  as  averaged  coupling 

U  tl 

functions  between  the  field  point  at  £  and  the  source  point  at  £'.  Area  1 
in  figure  5  corresponds  to  the  situation  in  which  both  the  field  point  and 
the  source  point  are  on  the  triangular  plate  S^,  while  area  2  corresponds  to 
that  in  which  they  are  both  on  the  rectangular  strip  S2*  For  areas  3  and  4 
the  two  points  are  on  different  sections  of  the  conductor. 

Of  the  four  areas,  area  2  is  the  most  important,  being  the  most 
extensive.  By  equation  33  rQ  and  are  identical  in  area  2: 

rD(£,£')  -  rH(£,£')  -  W(£- £')  for  £,£'  >  d  (56) 

By  equations  14,  39  and  40,  one  has 
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M.  nil!  . imp*  i.  i 


Figure  5.  Division  of  the  First  Quadrant  of  the  £-£'  Plane  into  Four 
Areas  of  Integration. 
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w(e-e')  =  ^5- 


fa 


4a 


dy 


fa 


d?'  At 


-jk/<5-  C')2  +  (y-y')2 


/(£-  £')2+  (y-y’)2 

-jk/(C-  C')2  +  Cy-y')2  +  4h2 
/(C-  5’)2+  (y-y’)2  +  4h2 


(57) 


Although  introduced  initially  in  area  2,  this  function  W  can  certainly  be 
continued  to  all  the  other  three  areas  in  the  £-£'  plane. 

For  the  evaluation  of  Fourier  transforms  it  is  convenient  to  decompose 
fp  and  rH  as  follows: 


rD(5,5')  -  W(£  -  £')  +  Ud(C,5') 

rH(c,cf)  -  W(£-  £')  +  uh(5,c’) 


(58) 


By  equation  56  the  two  functions  Up  and  vanish  identically  in  area  2. 

In  the  other  three  areas  they  are  both  expressible  as  a  difference.  Each 
of  them  has  a  very  complicated  form  and  will  be  written  out  in  full  in  the 
appendix.  The  advantage  of  the  decompositions  in  equations  58  lies  in  the 
fact  that  the  dominant  part  W  is  relatively  simple  and  continuous  throughout 
the  entire  range  of  integration  in  the  £-£'  plane.  The  discontinuities 
in  Tp  and  are  left  behind  in  Up  and  U^. 

The  Fourier  transforms  F  and  F  of  T  and  T  have  the  following 

u  HD  H 

corresponding  decompositions: 


where 


wD(P.P') 

wH(p»P’) 


fd(p.p’) 

fh(p,p’) 


wD(p,p’)  +  uD(p,p') 
wH(p,p’)  +  uH(p,p’) 

CD 

d£'W(£  -  £')sin  p£  sin  p'£' 

) 

d£  f  d£'W(£  -  £')cos  p£  cos  p'£' 

0 


(59) 


(60) 
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and 


u^p.p’) 


.00 

d£'U  (£,£’)sin  p£  sin  p'£' 

0 


^(p.p') 


d£'U  (£,£')cos  pC  cos  p’C' 


(61) 


It  will  be  shown  below  that  both  Wp(p,p')  and  w^(p,p')  contain  a  part 
proportional  to  6(p-p'). 

The  function  W(£ -£■'),  as  defined  in  equation  57,  can  be  converted 
to  a  form  more  suitable  for  Fourier  transformation.  This  is  achieved  by 
applying  the  Fourier  integral  representation  of  the  free-space  Green's 
function  given  in  equation  15.  The  integrand  in  equation  57  is  essentially 
made  up  of  free-space  Green's  functions.  Expressing  the  integrand  in 
equation  57  as  a  Fourier  integral,  one  can  immediately  perform  the  y- 
and  y' -integrations.  The  result  is 


W(C-  V) 


i-«-J2Th 


2  2  2  ,2 
a  +  3  +  Y  -  k 


(62) 


Furthermore,  one  can  carry  out  the  a- integration  and  obtain 

,2  ,  -j  2yh  ,J1  72  2*, 

1  -  e  J  -  j  /k  -  3  -  Y  I 


W(£-  Z') 


87r2j  •>_, 


d3 


) 


5-5’ 


00  —00 


X2~2  2 

/k  -  3  -  Y 


(63) 


where  the  square  root  is  defined  to  make 

/k2  -  32  -  Y2  +  k 


(64) 


as  k  goes  to  «.  Equation  63  offers  a  more  convenient  representation  of  the 
function  W(£-£')  than  equation  57,  since  the  £-  and  £' -dependences  are  all 
contained  in  a  simple  exponential  function.  The  Fourier  transform  of  W  can 
be  obtained  by  first  calculating  the  Fourier  transform  of  this  exponential 
function. 
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Consider  firsc  the  double  Fourier  sine  transform.  One  can  show  by 
straightforward  integration  that 


r«  e-l/k2-62-Y2|5-5'l8ll,  p£ 


-21  t^k.2  -  g2  - 

~2~2  '2 - I  sin 

k  -  B  -  Y  -  P 


.  2  „2  2  2 

k  -  g  -  Y  -  p 


-jA2  -  B2  -  y2  V 


(65) 


In  deriving  this  expression  one  must  keep  in  mind  that  k  has  a  small  negative 
imaginary  part,  as  Indicated  in  equation  16,  so  that  the  integrand  vanishes 
at  infinity.  Note  especially  that  the  right-hand  side  contains  a  term 
proportional  to  sin  p£'.  Multiplying  both  sides  of  equation  65  by  sin  p'C' 
and  integrating  over  £'  from  0  to  °°,  one  obtains  explicitly  the  following 
double  Fourier  sine  transform: 


f  d£  f 45  -62- v2ls- E'l 

'n  Jo 


sin  p£  sin  p’51 


FT ~2  2 

* j  2 — ~2 - T* — 2  5  <p  “  p ' )  +  n - 2 

k  -  6  -  y  -P  '  ’  - 


P2_ 


j ,  2  2  2  2  \/  2  „2  2  ,2\ 

\  k  -  6  -  Y  -P  A  k  -6  -y  -  p  J 


(66) 


The  6-function  comes  about  as  a  consequence  of  equations  47. 

The  double  Fourier  cosine  transform  can  be  evaluated  in  a  similar  way. 
The  result  is 


d5 


JO 


dC'e"^  "  3  "  Y  I5  “  Icos  p£  cos  p'C’ 


/I  12  2 

-rrj  -2  ~2  ^  2  5(P_  P'^  +  7~J 


,2  02  2 

k  -  S  -  y 


k  -  6  -  Y  -  P 


1 ,2  .2  2  2U.2  .2  2  ,2 

\k  -  S  -Y  -  P  k  -  g  -  y  -  P  - 


(67) 


Note  that  the  coefficient  of  the  6-function  in  equation  67  is  the  same  as 
that  in  equation  66. 
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Having  isolated  the  6-function  in  the  Fourier  transforms  in  equations 
66  and  67,  one  is  in  the  position  to  write  down  explicit  expressions  for  the 
functions  p(p)  and  A(p,p')  defined  in  equations  53  and  55.  The  result  for 
p(p)  is 


P<P)  “  — 7 
2  IT 


j.  /  sin  Ba  \2  1  -  e 


d6 


(68) 


The  y-integration  can  be  further  carried  out,  and  one  obtains 


2X  _  fi-j2^2-p2-g2h 
/k2  -  p2  -  B2 


(69) 


The  result  for  A(p.p')  is 


A(p.p’)  - 


2/j 


d6 


dY 


sin  Ba  \2  1  -  e'j2yh  k2fk2-82-Y2)  -  p2p'2 


Sa 


Y  1  -  e -  k  ik~-6~-Y~/  -  p-p' 

/  2  .2  2  (.2  .2  2  2  \(  ,2  2  2  ,2 

/k  -  6  -  y  '  k  -6  -y  ~P  /'  k  -6  -y  -P 


+  k  uH(p,p')  -  pp'uD(p,p') 


(70) 


If  one  wishes,  one  can  write  out  explicitly  the  two  functions  u^  and  u^ 
with  the  help  of  equations  61  and  the  appendix.  They  are  made  up  of  a  set 
of  complicated  four-dimensional  integrals. 


SECTION  VIII 


MODAL  STRUCTURE  OF  EXCITATION 


In  this  section  the  properties  of  the  integral  equation  54  will  be 
examined.  The  findings  will  not  only  provide  valuable  qualitative  informa¬ 
tion  on  the  solution,  but  will  also  suggest  the  appropriate  method  for 
constructing  an  approximate  solution. 

A  very  common  technique  for  solving  an  integral  equation  of  the  type 
of  equation  54  is  iteration.  Using  this  technique  one  can  immediately  write 
down  the  following  first-iterated  solution: 


A(p) 


M(P) _ 

'  k2  -  p2  ) p(p) 


(71) 


where 


M(p) 


2^eoVo 

2 


A(p.p') 

(  k2-p’2  )  p(p’)  . 


(72) 


Although  this  approximate  solution  is  not  expected  to  be  very  accurate, 
and  will  not  be  used  in  the  quantitative  calculations  of  this  study,  it 
does  reveal  many  of  the  general  properties  of  the  true  solution. 

The  behavior  of  the  solution  in  equation  71  is  dominated  by  the  zeros 
of  the  denominator.  These  zeros  give  rise  to  infinities  in  A(p) .  Each  of 
the  infinities  corresponds  to  a  discrete  mode  of  electromagnetic  propagation 
along  the  simulator  structure  in  figure  2.  Besides  these  discrete  modes, 
which  are  infinite  in  number,  the  solution  also  exhibits  a  spectrum  of 
continuous  modes,  as  will  be  seen  below. 

The  first  two  zeros  of  the  denominator  are  easily  found  to  be  at 
p  *  t  k.  They  correspond  to  the  TEM  mode  of  propagation.  All  the  other 
zeros  are  contained  in  the  function  p(p).  Inasmuch  as  p(p)  is  presented 
in  the  form  of  a  definite  integral  in  equation  69,  these  zeros  are  not  as 
easy  to  locate  accurately.  For  the  purpose  of  a  qualitative  investigation, 
however,  one  can  use  an  approximate  evaluation  of  the  integral.  When  the 


i 

f 


Is 


half-width  a  of  the  rectangular  strip  can  be  considered  large  relative  to 

2 

the  height  h,  the  factor  (sin  3a/8a)  in  the  integral  in  equation  69,  as 
compared  with  the  remainder,  is  sharply  peaked  at  £*0.  As  an  approximation, 
one  can  replace  the  remainder  by  its  value  at  6  =  0,  and  obtain 


p(p) 


-j/k2"^ 

~  a  J 


p  h  sin1 


(/ k2  -_e3  h  )  1 


/k2  - 


dB 


^  sin  Ba 


6a 


): 


(73) 


*  Thezeros  of  p(p)  in  equation  73  are  all  contained  in  the  factor 
sin  (  /k2  -  p2  h  )  p^",  which  is  an  entire  function  of  p.  It  has  an 

infinite  number  of  simple  zeros  located  at 


n  =  1,2,3, .. . 


(74) 


Note  that  because  of  the  square  root  in  the  denominator  of  the  factor  the 
case  n  -  0  is  not  included.  These  values  of  p  correspond  to  the  propagation 
constants  of  the  higher-order  modes.  They  are  seen  to  be  either  purely 
real  or  purely  imaginary,  implying  thereby  that  some  higher-order  modes 
are  unattenuated.  However,  if  one  had  used  equation  69  instead  of 
the  approximate  equation  73  to  compute  the  zeros,  one  would  have  obtained 
complex  propagation  constants.  One  would  have  concluded  that  all  higher- 
order  modes  are  more  or  less  damped. 

In  addition  to  the  zeros,  p(p)  has  branch  points  at  p  =  ±  k,  as  the 
exponential  factor  in  equation  73  shows.  The  branch  cuts  of  p(p)  in  the 
complex  p-plane  correspond  to  a  continuous  spectrum  of  propagation  constants. 
The  existence  of  continuous  modes  is  typical  of  open  waveguide  structures. 

Thus,  based  on  the  analytical  properties  of  the  denominator  in 
equation  71,  one  concludes  that  the  excitation  on  the  simulator  can  be 
decomposed  into  certain  discrete  and  continuous  modes.  One  can  show  that 
these  modes  are  in  fact  the  propagation  modes  of  an  infinite,  parallel-strip 
waveguide.  The  propagation  constants  of  the  lowest  modes  of  this  waveguide 
have  been  computed  accurately  in  reference  4.  The  simple,  approximate 
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formula  for  the  higher-order  propagation  constants  in  equation  74  compares 
very  well  with  the  accurate  values. 

From  these  findings  it  is  clear  that  the  free  term  on  the  left-hand  side 
of  the  integral  equation  54  describes  the  properties  of  a  parallel-strip 
waveguide.  The  role  of  the  integral  with  kernel  A(p,p')  is  to  bring  into 
consideration  other  aspects  of  the  simulator  geometry.  According  to 
equation  70,  A  is  made  up  of  two  parts.  The  first  part  is  written  out 
as  a  two-dimensional  definite  integral,  and  does  not  involve  the  linear 
dimension  d  of  the  triangular  plates.  Its  occurrence  in  the  kernel  is  a 
result  of  the  fact  that  the  simulator  is  a  semi-infinite,  rather  than  an 
infinite,  waveguide.  The  second  part  of  A  in  equation  70  involves  the 
parameter  d,  and  further  elaborates  the  simulator  geometry.  It  describes 
the  fact  that  the  simulator  has  a  discontinuity  in  the  form  of  a  junction 
between  a  section  of  a  conical  waveguide  and  a  parallel-strip  waveguide. 

In  equation  71  these  effects  of  A  enter  in  the  numerator  M(p) .  They  affect 
the  amplitudes  and  phases  of  the  modes  which  make  up  the  complete  solution. 
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SECTION  IX 


SELF-CONSISTENT  SOLUTION 

The  recognition  that  the  solution  is  made  up  of  modes  with  well-defined 
propagation  constants  suggests  a  self-consistent  method  for  solving  the 
integral  equation  54.  In  this  method  one  first  expresses  the  total  current 
flowing  on  the  simulator  as  a  linear  combination  of  the  propagation  modes 
of  the  parallel-strip  waveguide.  Each  mode  is  represented  by  a  simple 
exponential  function  of  5  involving  a  definite  propagation  constant.  Such 
an  expansion  contains  a  number  of  unknown  expansion  coefficients.  One  then 
evaluates  the  Fourier  cosine  transform  of  this  expression  and  substitutes 
it  into  the  integral  in  equation  54  as  a  trial  solution  in  place  of  the 
unknown  function  A(p').  The  resulting  equation  can  be  solved  for  A(p) , 
and  the  solution  is  of  the  form  in  equation  71.  As  one  has  already  seen, 
this  solution  itself  is  again  a  linear  combination  of  the  propagation  modes 
of  the  parallel-strip  waveguide.  Consequently,  the  expansion  coefficients 
can  be  determined  self-consistently . 

A  self-consistent  calculation  can  be  performed  with  an  arbitrary  degree 
of  elaborateness,  depending  on  the  number  of  expansion  coefficients  entered 
in  the  trial  solution.  In  what  follows,  the  integral  equation  54  will  be 
solved  using  the  simplest  self-consistent  scheme.  That  is  to  say,  the 
trial  solution  employed  will  involve  only  the  TEM  mode  of  propagation. 

In  view  of  the  geometry  of  the  simulator  and  the  manner  in  which  it  is 
excited,  the  current  will  be  predominantly  in  the  TEM  mode.  Hence,  a  TEM 
trial  solution  is  expected  to  be  satisfactory.  One  must  emphasize,  however, 
that  this  in  no  way  means  that  the  resulting  solution  exhibits  only  a  TEM 
mode.  It  means  simply  that  only  the  TEM  component  is  determined  self- 
consistently. 

The  current  associated  with  the  TEM  mode  is  given  by 

IT  (C)  =  I  e~iU  (75) 

TEM  o 
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where  the  amplitude  I  is  to  be  determined  self-consistently .  The  Fourier 
cosine  transform  of  the  TEM  current  can  be  easily  evaluated: 


dC  ITEM(C)cos  pC 


2jk 

IT 


(76) 


on  account  of  equation  16.  Substituting  the  right-hand  side  of  equation  76 
into  the  integral  of  equation  54  as  a  trial  solution  for  A(p'),  one  obtains 
the  following  approximate  solution  of  the  integral  equation: 


where 


A(p) 


N(P) 

(  k2  -  p2 )p(p) 


N(p) 


2jaie  V 
o  o 

2 

7T 


A(p>p') 

,2  ,  2 

n  1  -  k 


(77) 


(78) 


The  total  current  is  given  by 


1(5) 


.00 

0 


cos  pC 


(79) 


This  expression  provides  an  approximate  solution  of  the  problem  up  to  an 

unknown  constant  I  .  One  then  determines  this  constant  self-consistently 
o 

by  separating  out  the  TEM  part  of  1(C)  and  substituting  it  into  the  left- 
hand  side  of  equation  75. 

It  can  be  shown  that  both  N(p)  and  p(p)  are  even  functions  of  p,  so 
that  equation  79  can  be  rewritten  as 


KO 


dp 


e 


-jpC 


(80) 


With.  1(C)  in  this  form  one  can  evaluate  the  TEM  part  by  the  method  of  residues 
The  value  of  the  integral  is  dependent  on  the  singularities  of  the  integrand 


in  the  complex  p-plane.  As  has  been  discussed  in  section  VIII,  the 
denominator  of  equation  80  gives  rise  to  simple  poles  and  branch  cuts. 

They  correspond  to  the  discrete  and  continuous  modes  of  the  parallel-strip 
waveguide.  These  singularities  are  shown  in  figure  6.  The  singularities 
of  the  numerator  N(p),  if  any,  need  not  be  considered  here. 

The  coordinate  £  is  always  positive.  For  C  sufficiently  large,  the 
path  of  integration  in  equation  80  can  be  closed  by  a  large  semicircular 
arc  in  the  lower  half  of  the  complex  p-plane.  Then  the  TEM  part  of  1(0 
is  given  by  the  residue  at  the  pole  p  =  k: 


ITEM(5) 


j  7tN  (k)  -jk£ 
2kp (k)  6 


(81) 


Substituting  aquation  81  into  equation  75  and  using  equation  78  one  can 
solve  for  I  ,  and  obtain 


where 


(82) 


(83) 


With  the  value  of  Iq  determined,  equation  80  becomes  an  explicit,  approximate 
formula  for  the  total  current  flowing  on  the  simulator. 

It  is  possible  to  improve  systematically  the  accuracy  of  the  solution 
by  including  more  and  more  higher-order  modes  in  the  trial  solution. 

The  physical  content  of  the  self-consistent  solution  in  equation  82 
can  be  more  easily  appreciated  if  one  rewrites  it  in  the  following  form: 


I 

o 


V 

_ o 

z1  +  z2  +  z 


3 


•  (84) 
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where  Z^,  Z 2  and  Z3  are  three  quantities  with  the  dimension  of  an  impedance. 
They  are  defined  by 


Z1  m  irp(k)Zo 


—CO 


dg 


( 


sin  ga  \ 

Sa  ) 


2 


(85) 


and 


(86) 


(87) 


Z2  and  Z^  represent  the  integral  of  A  in  equation  82.  Z2  comes  from  that 
part  of  A  that  is  expressed  as  a  two-dimensional  integral  in  equation  70, 
while  Z3  comes  from  the  remainder.  These  impedances  are  functions  of  the 
three  dimensionless  parameters  h/a,  kh  and  d/h  of  the  problem.  It  can  be 
shown  that  Z^  depends  only  on  h/a,  Z2  depends  on  h/a  and  kh,  and  Z^  depends 
on  h/a,  kh  and  d/h: 

h-hii) 

h-h  (!•“*)  <88) 

Z3  ’  Z3  (  a  ’  kh’  h  ) 


Equation  84  interprets  the  TEM  current  Iq  as  equivalent  to  that  resulting 
from  the  application  of  a  voltage  VQ  across  the  terminals  of  three  impedances 
Z^,  Z2  and  Z^  connected  in  series.  The  physical  origin  of  these  impedances 
can  be  traced  among  the  findings  of  section  VIII.  Z^  is  the  characteristic 
impedance  of  an  infinite  parallel-plate  transmission  line.  Its  expression 
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in  equation  35  is  of  course  not  exact,  but  obtained  under  the  approximation 
of  averaging  over  the  transverse  direction.  Z2  is  a  correction  to  Z^,  and 
takes  account  of  the  fact  that  the  simulator  model  considered  in  this  study 
has  a  semi-inf inite,  rather  than  infinite,  structure.  is  an  equivalent 
impedance  for  the  junction  between  the  two  different  sections  of  the 
simulator,  and  embodies  the  effects  of  the  junction  discontinuity.  It 
alone  depends  on  the  dimension  d  of  the  conical  transmission  line  section. 

The  integral  in  equation  85  can  be  evaluated  analytically.  One  obtains 


and,  with  the  help  of  equations  61,  equation  87  can  be  rewritten  as 


Z  -  2jkZ  d? 

j  o 


0  '0 


,00  r 
ds’ 


cos  k£  Uh(£,£’)  -  j  sin  k£ 


VE*E’>] 


where  the  functions  0  and  Un  are  provided  in  the  appendix. 


SECTION  X 


NUMERICAL  RESULTS 


This  section  presents  the  results  of  numerical  calculations  on  the 
three  impedances  Z^,  Z^  and  In  these  calculations  the  following  sets 

of  values  of  the  three  parameters  in  the  problem  are  considered: 


-  -  1,  2,  3 
a 


kh  =  0.1,  1,  10 


(91) 


^  »  2,  4,  8,  12,  16,  20 


with 


L  *  d  cos  0 


(92) 


L  can  be  regarded  as  the  length  of  the  conical  transmission  line  section. 
The  ratio  d/h  is  related  to  L/h  by  the  simple  equation 


t-AIf 


+  i 


(93) 


For  a  typical  value  of  50  meters  for  h,  the  three  chosen  values  of  kh  in 
equation  91  correspond  approximately  to  frequencies  of  100  kHz,  1  MHz  and 
10  MHz,  respectively. 

The  impedance  Z^  is  purely  real,  and  depends  only  on  the  ratio  h/a. 

It  can  be  computed  easily  from  equation  89.  The  results  are  shown  in 
table  1.  As  has  been  pointed  out  earlier  in  section  IX,  Z^  is  in  fact 
the  characteristic  impedance  of  an  infinite,  parallel-plate  transmission 
line.  It  has  been  investigated  extensively  by  the  method  of  conformal 
mapping  in  references  5  and  6.  The  numerical  values  of  Z^  obtained  in 
reference  6  are  quoted  in  table  1  for  comparison.  A  glance  at  the  table 
reveals  that  the  two  sets  of  impedance  values  calculated  by  the  two  totally 


Table 


different  methods  are  in  uniform,  good  agreement  with  each  other.  The 
difference  is  only  about  5%.  The  agreement  is  also  found  to  be  equally 
good  for  values  of  h/a  other  than  the  three  listed. 

The  numerical  results  for  the  impedance  are  shown  in  table  2. 

It  is  seen  that  the  real  part  of  Z ^  is  negative.  But  the  combination 
Z^+Z£,  which  can  be  regarded  as  the  characteristic  impedance  of  a  semi¬ 
infinite  parallel-plate  transmission  line,  is  passive  as  it  should. 

Table  3  shows  the  numerical  results  for  the  impedance  Zy  They  are 
obtained  mostly  for  kh  =  0.1.  This  is  because,  at  higher  values  of  kh,  the 
integrand  in  equation  90  contains  a  large  number  of  oscillations.  The 
numerical  calculation  then  becomes  very  time  consuming.  The  table  does 
include,  however,  two  cases  at  kh *  1  and  2.5,  and  at  values  of  h/a  and  L/h 
that  are  comparable  to  the  dimensions  of  an  existing  simulator. 

At  kh  «*  0.1,  table  3  shows  that  Z^  is  mainly  imaginary.  Its  magnitude 
decreases  with  increasing  L/h  for  all  three  values  of  h/a.  Based  on  these 
data  one  can  say  that  the  transition  across  the  junction  is  smoother  the 
larger  the  ratio  L/h  becomes. 
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Table  3 


VARIATION  OF  IMPEDANCE  Z 3  WITH 

PARAMETERS  kh,  h/a  AND  L/h  (Z  -  4T7e  *  377  OHMS) 

0  0  0 


SECTION  XI 


DISCUSSIONS 


In  this  study  one  has  succeeded  in  obtaining  an  approximate  analytical 
formula  for  the  total  electric  current  flowing  on  a  parallel-plate  EMP  simulator. 
This  formula  enables  one  to  examine  both  qualitatively  and  quantitatively  the 
nature  of  the  excitation  on  the  simulator.  Contacts  with  previous  work,  on 
different  aspects  of  the  simulator  have  been  established.  The  formula,  however, 
cannot  yield  directly  the  electromagnetic  fields  in  the  simulator's  test 
volume.  In  order  to  determine  the  field  components  one  will  have  to  make 
further  approximations. 

One  approximation  consists  of  synthesizing  the  fields  in  the  test 
volume  from  the  field  patterns  of  the  discrete  modes  of  an  infinite,  parallel- 
plate  waveguide.  These  field  patterns  have  been  calculated  in  reference  4. 

The  fields  in  the  test  volume  can  be  expanded  in  terms  of  these  patterns. 

It  remains  to  determine  the  expansion  coefficients.  As  an  approximation 
the  coefficients  can  be  taken  as  identical  to  the  expansion  coefficients  of 
the  total  current  in  terms  of  the  modal  currents.  That  is  to  say,  they  are 
given  by  the  residues  of  the  integrand  in  equation  80  at  the  poles  in  the 
complex  p-plane  corresponding  to  the  discrete  waveguide  modes. 

Another  method  consists  of  constructing  approximately  an  expression  for 
the  surface  current  density  from  that  for  the  total  current  in  equation  80. 

The  variation  of  the  surface  current  density  on  the  simulator  In  the  transverse 
direction  is  largely  determined  by  the  edge  condition.  One  can  distribute  the 
total  longitudinal  current  across  the  breadth  of  the  simulator  according  to 
some  simple  analytical  function  which  satisfies  the  edge  condition.  The 
electromagnetic  fields  can  then  be  obtained  from  the  surface  current  density 
by  integration. 
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APPENDIX 


FUNCTIONS  Ud(5.S'>  AND  UH(C,C) 

The  functions  U_.(5,£')  and  U„(C,5')  first  occur  in  equations  58.  They 
D  £1 

can  be  calculated  from  the  relations 


Ud(5,5')  =  -  W(S-e') 

uh(5,s')  =  rH(c,e')  -  w(5-c') 


(Al) 


rQ  and  Fh  are  certain  averages  of  the  Green's  functions  GD  and  G^  defined  in 
equations  14,  31,  32  and  33,  while  W  is  given  in  equations  56  and  57. 

U..  and  U„  have  different  forms  in  the  four  different  areas  of  the 
D  n 

plane  shown  in  figure  5: 


1.  In  area  1  (0  <  5,$’  <  d)  one  has 

e-jk/(c-c’)2  +  (y-y’)2 
/(C-C’)2  +  (y-y')2 

e-jk^2+C,2-2?C’cos  29  +  (y-y’)2 
42  +  C2- 2«'cos  28  +  (y-y')2 

e-jk/(5-C')2  +  (y-y')2 
/u-c')2  +  (y-y’)2 

-jk/?2  + ?'2  -  255'cos  26  +  (y-y')2 
-  cos  26  —  - - -  - - . 

/c2  +  C'2-  2C5’cos  29  +  (y-y1)4- 


GH(5,y,C',y') 


GD(5,y,5',y')  =  4^ 


Then 


(A2) 
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4.  In  area  4  (d  <  £  <  ®, 


0  <  5'  <  d)  one  has 


GD(C,y,Oy’) 


4r 


-jk/(5  -  d  +  d  cos  0  -  C'cos9  )2  +  (y  -  y ' )  2  +  (h  -  C '  sin  9) ‘ 


(_  /(£  -  d  +  d  cos  9  -  5' cos  9) 2  +  (y-y')2+(h-  5  'sin  9) 2 


-jk/(5  -  d  +  d  cos  0  -  C' cos  3) 2  +  (y  -  y ' ) 2  +  (h+  5 ' sin  9)‘ 


/  2  9  ,  2 

/(5  -  d  +  d  cos  9  -  5'cos  8)  +  (y  -  y ' )  “  +  (h  +  C  sin0  ) 


GH(?,y,?’ ,y')  -  cos  9  GD(?,y,5' ,y') 
Then 

d 


(A73 


U  (5,0  =  — J— 
u  /Or' 


U„(5,0  -  — ~ 

H  4a  V 


fa 

fa5 '  /  d 

dy'GD(5,y,5’,y')  - 

dy 

| 

-a  J 

-a  5*  /d 

(■a 

■a5'/d 

dy 

dy ' GH ( 5 » y , ? ' , y ' )  - 

-a 

-aC'/d 

(A8: 


Note  that  the  forms  of  U^U.O  and  UH(C,0  in  area  4  are  obtainable  from 
those  in  area  3  by  interchanging  the  variables  5  and  5'. 
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